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Abstract 

Using temperature dependent Raman spectroscopy we address the question of how the transition 
from bulk to few atomic layers affects the charge density wave (CDW) phases in lT-TaS 2 * We 
find that for crystals with thickness larger than lOnm the transition temperatures between 
the different phases as well as the hysteresis that occurs in the thermal cycle correspond to the 
ones expected for a bulk sample. However, when the crystals become thinner than « lOnm, the 
commensurate CDW phase is suppressed down to the experimentally accessible temperatures. In 
addition, the nearly commensurate CDW phase is diminished below ~ 4 nm. These findings suggest 
that the interlayer coupling plays a significant role in determining the properties of CDW systems 
consisting of a few unit cells in the vertical direction. 
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Transition metal dichalcogenides are typically layered materials, TX 2 (T= transition 
metal, X=chalcogen: S, Se,Te) built out of stacked three-layers units (X-T-X) coupled by 
weak van der Waals forces. Their bulk electrical properties span a wide range: from insulat¬ 
ing (HfS 2 ), semiconducting (M 0 S 2 ) to metallic showing superconductivity or charge density 
order (e.g. NbS 2 , TaS 2 , TiSe 2 ). The competition of various effects such as dimensional¬ 
ity, electron-phonon coupling, electron-electron interactions, and disorder gives rise to the 
wealth of phases in these materials [1], The anisotropy of the crystal and the weak layer 
bonding are typically understood as responsible for the quasi 2D properties of these mate¬ 
rials. However, attesting to the important role of the interlayer coupling, it was shown that 
when thinned down from bulk to only a few atomic layers some of these material properties 
change. For example, M 0 S 2 is in bulk an indirect-band semiconductor while one layer of 
the same material has a direct band gap [2]. Experimentally, exploring the questions of how 
thickness can be used to tune material properties is possible due to advances in fabrication 
techniques following the demonstration of micromechanical cleavage of atomically thin ma¬ 
terials [3]. In particular, the precise nature of the driving mechanisms as well as the nature 
of the charge density wave state in atomically thin layered materials is yet to be understood. 

The IT polymorph of TaS 2 , lT-TaS 2 , has one of the richest phase diagrams: it is metallic 
at higher temperatures, has four temperature-dependent charge density wave phases with 
different structures; under pressure and doping it becomes superconducting [4] and it was 
suggested to show Mott insulator behavior [5]. In this letter we use variable temperature Ra¬ 
man spectroscopy measurements to show that by reducing the thickness of lT-TaS 2 between 
~ 5 nm — 10 nm, the transition to a commensurate charge density wave is suppressed down 
to the experimentally accessed temperatures. Furthermore, the Raman spectrum suggests 
that for the crystals with thickness below ~ 4nm the nearly commensurate phase is absent. 

The structure of the lT-TaS 2 crystals is such that tantalum is octahedrally coordinated 
to the sulfur. In the c-axis, one unit cell is represented by a three-layer S-Ta-S as sketched 
in Figure 1(a). The different CDW phases that arise upon cooling and warming respectively 
are illustrated in Figure 1(b). Cooling down from ~ 600/t the system develops an incom¬ 
mensurate CDW (ICCDW) state and further cooling to approximate 350K takes it into a 
nearly commensurate CDW (NCCDW) arrangement. The term incommensurate refers to 
the fact that the CDW wavevector is not an integer of the atomic lattice wavevector. In the 
NCCDW, commensurate domains coexist with incommensurate ones. The transition at an 
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FIG. 1. (a) Schematic representation of the crystal structure for lT-TaS 2 - (b) Illustration of the 
different phases in lT-TaS 2 as a function of temperature, (c) Ta atoms arranged in the Star of 
David pattern in the CCDW phase, (d) STM topographic images of a bulk crystal of lT-TaS 2 
(T=59K, V Bias = 20 mV, I tunneling = 2nA) (e) Temperature dependence of the resistance of the 
lT-TaS 2 device indicated in the inset for both the cooling (blue) and warming (red) cycles. 

even lower temperature, ~ 150/v, is to a commensurate CDW (CCDW). In this situation, 
the Ta atoms participating in the associated periodic lattice distortion (PLD) are arranged 
in groups of 13 atoms, such a configurations being typically referred to as the Star of David 
(Figure 1(c)). When warmed up, the system shows hysteretic behavior so that it persists 
into the commensurate state until ~ 215 K as illustrated in Figure 1(b) [6]. 

Flakes thicker than 20 nm were characterized using resistivity measurements [7-9]. The 
samples were obtained by mechanical exfoliation from a bulk crystal of TaS 2 onto the surface 
of a Si/SiO 2 wafer, similarly to graphene exfoliation. An optical microscope was used 
to identify flakes of different thicknesses and atomic force microscope measurements were 
performed to determine the thickness. Typically, we obtained flakes with heights ranging 
from 2 nm to over 150nm. The devices were fabricated using a Laser Pattern Generator, 
photoresist as a lithographic mask and the contacts were thermally evaporated Ti/Au. An 
optical micrograph image of such a four-terminal device is shown as an inset to Figure 
1(e). The resistance was recorded as a function of temperature as presented in Figure 1(e). 
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FIG. 2. Temperature dependent Raman spectroscopy of a bulk crystal lT-TaS 2 * (a) Raman spectra 
on a bulk sample at the indicated temperatures corresponding to different CDW phases (blue - 
CCDW, red - NCCDW, black - ICCDW). The red and blue curve were recorded during a cooling 
cycle, while the black one is on a warming cycle, (b) and (c) Raman intensity map as a function 
of temperature for a thick sample in a cooling and warming cycle respectively. The sharp changes 
in the spectra represent the transitions indicated below the maps, (d) Top: Schematic of thinning 
down a crystal from bulk to a few layers. Bottom: Summary of the transition temperatures 
measured for different samples as a function of sample thickness. The hatched area represent a 
regime where the behavior was found to be distinct from bulk. 


The arrows indicate both the cooling cycle (blue curve) and the heating cycle (red curve). 
In the cooling cycle, the two jumps in resistance correspond to the transitions: ICCDW- 
NCCDW ( « 350° ) and NCCDW-CCDW (« 150°) respectively. When warmed up the 
CDDW-NCCDW transition occurs at « 220°. 
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FIG. 3. (a) Raman intensity map taken during a cooling cycle for a thick flake (h ~ 160nm) 
indicated in (c) by number 1. (b) Raman intensity map taken during a cooling cycle for a thin 
flake (h « (8 ± 2) nm) indicated in (c) by number 2. (c) Top: Optical microscope image of the 
lT-TaS 2 flakes where the data in (a) and (b) were taken. Bottom: Atomic force microscope height 
profile talken along the dashed line indicated above, (d) Comparison at the indicated temperatures 
between the Raman spectra taken for areas 1 and 2. (e) Raman spectra for the low wavenumbers 
at indicated temperatures for the areas 1 and 2. 


Microscopic details of the CDW phases in lT-TaS 2 can be visualized by Scanning Tun¬ 
neling Microscopy (STM) [10]. Using a variable temperature ultrahigh vacuum STM the 
CDW on a freshly cleaved crystal was readily imaged (Figure 1(d)). When zooming in, 
both the atomic lattice as well as the CDW is visible as it is the case presented in Figure 
1(d) which was taken at 59K. We find the period of the CDW, in agreement with previous 
literature. We note that the surface becomes degraded after exposure to air for more than 30 
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minutes, therefore the best conditions for a tunneling experiment were achieved by cleaving 
the crystal in-situ in the UHV system. 

Both techniques described before complement each other in unveiling the bulk electronic 
and microscopic details of the CDW phases. However, they become challenging for measur¬ 
ing flakes with only a few atomic layers. Resistivity measurements are hampered by high 
contact resistance, while STM measurements on exfoliated thin flakes would also involve 
nanofabrication techniques that contaminate the surface. An experimental method that 
can give insight into the CDW phases without the need for further processing of the thin 
exfoliated flakes is Raman spectroscopy. In this work we study the changes in the Raman 
scattering that correspond to transitions of the crystals into CDW phases, therefore prob¬ 
ing the lattice vibrations for both the only few-layers-thick crystals as well as for the bulk. 
While bulk crystals were studied before [11-13], the experimental results presented here for 
atomically thin flakes have not been reported previously. 

For each selected area the Raman spectrum is recorded as a function of temperature. We 
use a Renishaw InVia system with the laser 514nm and typical power below lmW . The 
system is equipped with a temperature controlled stage that allowed varying the temperature 
of the sample from 80K to 400K. The laser spot typically focuses on a diameter of ~ 3/im. 
It is known that thin flakes could be damaged by heating while shining the laser light during 
the Raman measurement [14]. For the experiments here, we checked the dependence of the 
spectrum on laser power and we confirmed that heating the thin flakes can be neglected 
for the laser power that was used in the experiment. The same checks ensure that the 
temperature of the flake was not affected by shining the laser light on them. 

For a bulk system we measured the Raman spectra at three different temperatures T = 
403AT, T = 233 K and T = 93K respectively as shown in Figure 2(a) . They have distinct 
features and they correspond to the ICCDW, NCCDW and CCDW phase respectively. We 
note that for this data set, due to the experimental set-up, the intensity of Raman spectrum 
below ~100cm _1 is strongly diminished. In the NCCDW phase the most pronounced phonon 
modes are seen at 108 cm -1 , 124 cm -1 , 244 cm -1 , 265 cm -1 , 282 cm -1 , 300 cm -1 , 307 cm -1 , 
327 cm -1 , 385 cm -1 . At lower temperatures, for the CCDW phase, we observe that the 
peaks become sharper and new phonon modes appear. The appearance of new peaks can 
be understood as a consequence of dispersion curves folding into the zone center due to the 
reduced size of the Brillouin zone. In this state, because the unit cell of the CCDW phase 
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contains 39 atoms, one expects 117 phonon modes from which 19 A g + 19 E g are Raman 
active [13, 15]. In our data, the Raman spectrum shows peaks at: 61 cm -1 , 71 cm -1 , 79 
cm -1 , 100 cm -1 , 114 cm -1 , 128 cm -1 , 256 cm -1 , 267 cm -1 , 277 cm -1 , 284 cm -1 , 298 cm -1 , 
307 cm -1 , 322 cm -1 , 331 cm -1 , 367 cm -1 , 388 cm -1 , 397 cm -1 . The most pronounced ones 
of the above list are also indicated by arrows in Figure 2(a). 

To illustrate the changes in the Raman spectrum as the system goes through the differ¬ 
ent transitions, we plot Raman intensity maps as a function of temperature as presented in 
Figure 2(b) and (c). A sequence of spectra in a cooling cycle is presented in Figure 2(b). 
The change in the spectrum occurs at T ~ 150 K corresponding to the NCCDW-CCDW 
transition. In Figure 2(c), for the same area as in Figure 2(c), the temperature dependence 
of Raman spectrum was recorded upon warming up the system. The abrupt change cor¬ 
responding to the transition CCDW-NCCDW is now present at ~ 215AT. The hysteretic 
nature of this transition (Figure 1(e)) can, therefore, be confirmed using the Raman spec¬ 
troscopy experiment. We observe that the modes significantly affected by the transition 
NCCDW-CCDW are the ones around 100 cm -1 . 



Raman shift (cm' 1 ) Raman shift (cm' 1 ) 


FIG. 4. (a) Raman spectrum taken at different temperatures, at the position indicated in the optical 
image of the flake as area 3. (b) Raman spectra for the indicated flakes of different thickness. 

We used such temperature dependent Raman maps as a way to explore how the different 
CDW phases change when thinning down the system to a few atomic layers. A summary 
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of the experimental data is presented in Figure 2(d). While for flakes thicker than ~ lOnm 
the system behaves similarly to the bulk, in the regime below h = lOnrn the behavior of the 
flakes was found to be significantly altered as explained below. 

The inset of Figure 3(c) shows a mechanically exfoliated flake that has a thicker part 
(yellow color), h\ ~ 160nm and a thinner part (black color) /i 2~(8 ± 2 )nm together with a 
height profile obtained by atomic force microscopy. For the two different areas of a sample, 
Figures 3(a) and (b) present the Raman intensity map upon cooling. The thick flake of 
area 1 shows the abrupt change as the system transitions from NCCDW to the CCDW as 
explained earlier. In contrast, the thin flake of area 2 , shows no transition into a CCDW, 
but it remains in the NCCDW at the lowest temperatures that were reached experimentally. 
For further comparison, in Figure 3(d) we plot the Raman spectra obtained on the areas 1 
and 2 at the temperatures indicated as inset. Consistently, when we experimentally probe 
the wavenumbers below 100cm -1 as plotted in Figure 3(e) we find the same behavior: while 
the Raman spectrum of the thick flake changed below the transition temperature, the one 
for the thin flake remains essentially identical to the one at higher temperatures. 

We now focus on samples that were below ~ 4nm. Those are typically transparent when 
seen through an optical microscope as illustrated in the inset of Figure 4(a). We measured 
the Raman spectrum of theses atomically thin flakes both at room temperature as well as 
low temperature. In Figure 4(a) we plot the Raman spectrum taken for the flake indicated in 
the inset by area 3 at T = 273 K, 13 6K, 93 K. We find that this spectrum is distinct from the 
one that corresponds to the NCCDW from Figure 2 (a) at all three temperatures, indicative 
that that for such thin flakes the NCCDW phase is also absent. When we measure flakes 
of intermediate thickness such as area 4, as seen in Figure 4(b), the spectrum develops a 
stronger peak at ~ 244cm -1 , suggesting that as a function of thickness the systems undergoes 
novel CDW reconstructions. 

The suppression of the CCDW and NCCDW phases for the ultrathin samples is strongly 
connected to the periodicity of the vertical stacking corresponding to the two phases. For 
the NCCDW, Xray experiments have determined the presence of a 3 -layer period [16, 17]. 
In the case of the CCDW phase some reports propose a disordered stacking [18] while others 
suggest a period of 13 layers. Theoretically, various models explore what are the energetically 
favorable stacking arrangements in the different phases [19, 20]. The critical thickness found 
by our experiments to influence the formation of the CDW phases is therefore consistent 



with these periods. The diminished CCDW below ~ lOnm and NCCDW below ~ 4 nm can, 
therefore, be understood as an effect of the crystal thickness becoming comparable to the 
c-axis period of the respective phase. 

In conclusion, using temperature dependent Raman spectroscopy we address the question 
of how the transition from bulk to few atomic layers affects the charge density wave phases 
in lT-TaS 2 - We find two thickness where the CDW phases of the crystals are significantly 
altered. The NCCDW is suppressed below ~ 4nm and the CCDW below ~ lOnm respec¬ 
tively. The findings are consistent with the previously reported c-axis periods for the two 
phases. These experimental results attest to the importance of the interlayer coupling to 
determine the properties of ultrathin van der Waals materials and open up possibilities for 
tuning the phase diagram of materials by varying their thickness. 
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